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Abstract 

We have observed thermal gating, i.e. electrostatic gating induced by hot electrons. The effect 
occurs in a device consisting of two capacitively coupled quantum dots. The double dot system is 
coupled to a hot electron reservoir on one side (QDl), whilst the conductance of the second dot 
(QD2) is monitored. When a bias across QD2 is applied we observe a current which is strongly 
dependent on the temperature of the heat reservoir. This current can be either enhanced or 
suppressed, depending on the relative energetic alignment of the QD levels. Thus, the system can 
be used to control a charge current by hot electrons. 
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In recent years thermoelectrics, thermionics and thermal management in small-scaled 
devices have become important subjects both in basic and applied solid-state research BE]- 
In order to control the heat flow on the nanometer scale, a number of promising concepts 
have been proposed and partly realized, including solid-state thermal rectihers [3HZ], thermal 
transistors PE], and nano-refrigerators mm- Moreover, new concepts for highly efficient 
energy harvesting devices have recently been proposed which use a system of two quantum 
dots (QDs) as a central building block to convert heat into a directed current [T21IH] . The 
key feature of these systems is the capacitive inter-dot coupling [TS] which enables energy 
exchange between the QDs while particle transfer is blocked. Here we show how a system 
of two Coulomb coupled QDs acts as a thermal gate for charge currents. One of the dots 
(QDl) can exchange electrons with a hot reservoir only. The other dot (QD2) connects two 
reservoirs of equal but lower temperature. If a small potential difference is applied across 
QD2, we observe that the resulting current can be either enhanced or suppressed by variation 
of the temperature in the hot bath connected to QDl. An intuitive picture is given which 
explains the underlying mechanism. It is shown that this effect of thermal gating is in fact 
strongly related to correlations between fluctuations in the occupation number of both QDs. 

The device is processed by means of optical and e-beam lithography and subsequent 
metalization of Ti/Au-electrodes (gates) on a GaAs/AlGaAs heterostructure with a two- 
dimensional electron gas (2DEG) 94 nm below the surface (carrier density n = 2.4 x 10^^ 
cm“^, electron mobility /i = 0.69 x 10® cm^/Vs). The gate structure is shown in Fig. (a). 
The gates (black) are labeled 1 through 7, PI and P2. The electron reservoirs are denoted 
S, D and H. The QDs are dehned by gates 3 to 7, labeled QDl and QD2. QDl is in 
direct contact with reservoir H (red), which will serve as a hot electron reservoir. QD2 is 
connected to reservoirs S and D (both blue) which are at a lower temperature, representing 
the conductor circuit. The junctions are tuned into the tunneling regime by adjusting 
the voltages of the appropriate gates. The plunger-gate voltages PI and P2 control the 
chemical potentials and of the QDs such that the electron occupation numbers 
can be adjusted individually. The two QDs are separated by an electrostatic barrier (gate 
5), which suppresses the electron transfer between the dots. At the same time, QDl and 
QD2 interact electrostatically due to the small spatial separation. This means that the 
occupation number of each QD affects the chemical potential of the other QD, and thus 

M)] where N and M are the occupation numbers of QDl 
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and QD2, respectively mm- Reservoir H defined by gates 1, 2, 3 and 4 forms a channel 
of 20 /im length and 2 pm width. Opposite QDl, a constriction is created by gates 1 and 2 
which can be used as a voltage probe in the channel. Its conductance is set to G = lOe^/h 
by adjusting the gate voltages Vi and V 2 thus ensuring that no thermovoltage is created 
across this junction when the temperature in reservoir H is increased |18j . 

The sample is mounted in a top-loading dilution refrigerator with base temperature 
^base = 80 mK. For a conductance characterization of QD2 with all reservoirs at Tbase, 
an excitation voltage 14c = 5 pV (/ = 113 Hz) is applied between reservoirs S and D. The 
current is measured with a lock-in, using a current ampliher that connects D to a virtual 
ground potential. By varying Vpi and 14>2 one obtains the so-called stability diagram of the 
QD-system na. shown in Fig. (b) where the conductance G of QD2 is displayed in a gray 
scale as a function of the voltages Vpi and 14 > 2 - 

Along the horizontal axis 14>2, we observe two conductance resonances which identify those 
gate voltage conhgurations for which p^^^ is aligned with ps and po- They are separated by 
the Coulomb charging energy of QD2. Due to the mutual capacitive coupling the energetic 
position of p^^^ is affected by the energy of QDl. This leads to a continuous shift of the 
conductance resonances for larger 14>i towards smaller 14>2 [dashed, red lines in Fig. HKb)). 
When p(^^ aligns with pn, N changes by one [solid, red lines in Fig. (b)]. This causes 
discrete jumps for the conductance resonances of QD2, indicated by red arrows in Fig. (b). 
These jumps are a result of the capacitive inter dot coupling which leads to the transfer of 
the energy Eq'. -|- 1,M) = M) -f- Eq [151 [IZ]- Hence, the charge occupation 

numbers of both QDl and QD2 are stable only in the regions enclosed by solid and dashed 
lines in Fig.[^(b). The energy Eq can be calculated from the displacement of the conductance 
resonance along the l/p 2 -direction, AVc; indicated by yellow dotted lines in the hgure. Using 
the gate efficiency 02 = 0.032 obtained from d//dU characterization of QD2 yields Ec ~ 
90 pV. 

In order to subject the QD-system to a temperature difference, we make use of a current 
heating technique [TS]: An ac-current Jh = 150 nA with frequency / = 113 Hz is applied 
to the heating channel (reservoir H). Because of the strongly reduced electron-lattice in¬ 
teraction in GaAs/AlGaAs 2DEGs at low temperature, the energy is dissipated into the 
lattice only in the wide contact reservoirs. On a length scale of a few pm, however, electron- 
electron scattering dominates electron-phonon scattering, resulting in a thermalized hot 
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Fermi distribution of the electrons in the channel only. Based on QPC-thermometry [18] we 
estimate that for a current of Jh = 150 nA, Th increases by AT 100 mK. The ac-heating 
causes the temperature in the heat reservoir to oscillate at 2/ = 226 Hz between Tbase and 
^max = Tbase+AT. Tliis ensures that all temperature-driven effects also oscillate at frequency 
2/, enabling straight forward lock-in detection. Next, a dc-voltage source is connected to S 
which applies Vs,gnd while the current ampliher (input impedance i?imp = 2 kf2) connecting 
D to ground potential is read out by a lock-in ampliher detecting at 2/ = 226 Hz. This 
allows us to determine the change of the current in the drain contact AJd due to variation 
ofTn. 

With Vs,gnd ~ — 30/iV we obtain the data shown in Fig. (c). The lines delimiting 
the stability regions are indicated. Surrounding each stability region vertex we observe a 
four-leaf clover shaped structure that is composed of positive and negative current changes 
of up to ±8pA. The sign changes occur at the transitions from one quarter of a “clover 
leave” to the adjacent ones. Diagonally opposite regions exhibit identical sign. A closeup of 
a similar clover-leaf structure, obtained for slightly different values for Vpi and Vp 2 , is given 
in Fig. (a) for Vs,gnd = —100/iV. [For the measurements shown in Fig. the current 
ampliher is replaced by a resistor R = 100 kf2, the voltage drop across which is detected by 
the lock-in at 2/.] The corresponding conductance stability vertex is shown in Fig. (c). A 
direct comparison identihes the four parts of the clover-leaf pattern with diherent stability 
regions of the vertex: Sections 1 (A^+1, M) and 4 {N, M+1) produce a positive signal 
while for sections 2 (A-l-1, M+1) and 3 {N, M) negative A Id are observed. A single 
trace extracted from the color scale plot for constant Vpi = 381 mV (green, horizontal line) 
is shown in the top panel of Fig. (a). It exhibits a maximum and a minimum at Vpc 
corresponding to dot occupation (A-l-1, M) and (A^-l-l, M+1), respectively. In between, 
the signal changes approximately linearly with Vp 2 . Moving further away from the vertex 
causes the signal to decay. A trace extracted along the Vpi axis for constant Vp 2 = 521 mV 
(red, vertical line) behaves likewise (side panel). In a next step the dc-voltage applied to S 
is reversed, so that Vs,gnd = 100/iV. The result is given in Fig. (b). Clearly, the clover 
leaf pattern is reproduced, however, with all signs inverted. 

We now discuss qualitatively how we can understand this behavior. As is evident from 
Fig-i the sign of AJd does not change over a single stability section of the system [labeled 
1 - 4 in Fig. (a) and (c)]. Furthermore, the 2/-detection of the signal indicates that these 
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current signals are triggered by a temperature change in reservoir H. In the vertex region, 
the occupation numbers of both QDs can fluctuate while the occupation number becomes 
hxed when moving away from this region. It is thus apparent that the current changes which 
give rise to the clover-leaf structure originate from fluctuating occupation numbers of QDl 
and QD2. As an example, Fig. (a) shows the alignment of and with Vsd < 0 for 
section 1 with A^-|-l electrons on QDl and M electrons on QD2. Due to the ac-character of 
the heating current Th oscillates between the two values Th = Tbase and Th = Tmax > ^base- 
The hrst case is shown on the left side of Fig. (a): QDl is occupied with A^-|-l electrons, 
i.e. is below pn and therefore the electron number of QDl is hxed at A^-|-l. QD2 is 
occupied with M electrons and the chemical potential p*'^^(A^-|-l, M+1) which is required 
to add the electron lies outside the bias window Vsd- Thus, transport across QD2 

is blocked. When the temperature in reservoir H is increased such that Th = Tmax [right 
hand side in Fig. (a)], empty states are created below ph in this reservoir. This increases 
the charge huctuation rate on QDl. However, when due to these huctuations QDl relaxes 
to the A^-state, the energy required to add an electron to QD2 is reduced by Eq- The 
corresponding {N, M+1) is below ps and the current across QD2 increases [indicated 
by red arrows in Fig. (a)]. Since Th oscillates between Tbase and T^ax, this effectively 
leads to a temperature driven modulation of the conductance of QD2: If Th increases, the 
current across QD2 increases as well. For Th at a minimum, transport is blocked. The 
resulting current modulation at the drain contact is then detected by the lock-in ampliher 
as a positive signal. 

The QD levels for section 2, with occupations (A^+1, M+1), are depicted in Fig. |(b). 
Starting again with the condition Th = Tbase so that A^-|-l is hxed, we hnd that transport 
across QD2 is enabled because p*^^^(A^-|-l, M+1) is situated within I^d- However, charge 
huctuations on QDl, which increase with increasing Th [right side in Fig.|^(b)], tend to block 
transport across QD2: The corresponding M+1) is below pd and thus, electrons are 

trapped on QD2. The correlation between Th and Id is now inverted compared to section 1: 
a temperature increase tends to block transport while small Th increase Id- The sign of the 
signal in sections 3 and 4 can be understood in a similar manner: In these sections QDl is in 
the Wstate and huctuations lead to an occasional occupation with A^-|-l electrons. This is a 
reversal of the situation in sections 1 and 2 because an increase of Th now causes a decrease 
of Id when QD2 is in the M-state (section 3) while it enhances Id for QD2 exhibiting M+1 
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electrons (section 4). Within this picture, the explanation of the observations for a reversed 
voltage bias [I^d = 100 pV, cf. Fig. |^(b)] is also straightforward: Because a sign change of 
the bias voltage reverses the dc-current through QD2 this leads to an overall reversal of the 
observed signal. 

We have performed simple model calculations to substantiate the qualitative discussion 
presented above. Assuming sequential transport across QD2, the current Id can be related to 
the applied difference in electro-chemical potential Vsd = P- 5 —by considering Fermi-Dirac 
occupation statistics f{n^'^\Tj) = 1/(1 -|- exp(p*^^^ — Hj/k^Tj)), j = S,D in the source and 
the drain contact and a single resonant QD level which is located at p*^^^ = —Ec/2. For 
P 5 < Pd we can then write /^(p^^^) oc f{—Ec/2,Ts) x (1 — /(— i?c/2, Td)). The current 
across QD2 when QDl hosts A^-|-l electrons can be treated likewise, with = +Ec/2. The 
total current Id through QD2 is now the sum of and 7^, weighted with the appropriate 
probabilities of QDl hosting N or 77-1-1 electrons. Thus, J£)(p^^\ oc 

Figure (a) compares Id as a function of p^^^ and p*^^^ for Th = Tbase (left) and 
Tr = Tbase + AT = Tmax witli AT = 100 mK (right) while Ts,d = Tbase = 230 mK, A/i = 
100 pV and Eq = 90 pV. As expected, the results strongly resemble the conductance stabil¬ 
ity diagram in the vertex region. However, major differences for different Th are not directly 
obvious. In order to model our experiment we subtract the calculated data sets in Fig. (a) 
from each other such that we obtain AJd = 7D(Tmax) — 7 d(T base), which corresponds to 
the change in current through QD2 due to a change of Th by AT. The result is given in 
Fig. 1^ (b). Evidently, the clover-leaf pattern found in the experiments, is reproduced nicely. 

We point out that a similar four-leafed clover pattern has been observed in connection 
with Coulomb coupled double QDs previously: McClure et al. [I6] have reported on exper¬ 
iments addressing the cross-correlation of shot noise in such a system. There, the authors 
observed regions of positive and negative correlation at the stability region vertex arranged 
in a cloverleaf shaped pattern similar to the one discussed here. The underlying mechanism 
is actually closely related to the one active in our experiments: Negatively correlated shot 
noise indicates that charge fluctuations of one QD tend to suppress fluctuations on the other 
one (and vice versa). Correspondingly, those are the conhgurations for which we observe 
a reduced current through QD2 if the temperature in reservoir H is increased. A positive 
correlation implies that occupation fluctuations tend to occur simultaneously on both dots. 
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Thus, we observe an enhancement of the current through QD2 with temperature in those 
regions. 

In order to estimate the gating range of our device we analyze the data shown in Figs. (a) 
and|^ (c): Using the G of QD2 at those conhgurations for which a maximal A/d = 18 pA is 
observed {G = 0.09 e^/h) we calculate the drain current for Vsd = 100/iV and Th = Tbase, 
which gives Id = 360 pA. Relating this current to A/d then yields a gating amplitude of 5%. 
Although this ratio is rather small, it can be strongly enhanced by tuning the parameters 
Eq and /cbAT. 

Finally, we note that the thermal gating effect presented here could be used, e.g., to mon¬ 
itor carrier heating in quantum circuits. Furthermore, it could be utilized to also manipulate 
heat flow across QD2: Since the thermal conductance k of a QD as a function of /i usually 
follows the Wiedemann-Franz rule and thus has a similar line shape as the conductance [T9] , 
the mechanism presented here would allow gating of heat currents to be accomplished, thus 
suggesting a route to realizing a QD-based all thermal transistor. 

The authors gratefully acknowledge discussions with B. Sothmann, R. Sanchez and M. 
Biittiker. This work has been hnancially supported by the DFG (SPP1386). 
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FIG. 1. (a) Schematic design of the gate structure (black). Gates are labeled with numbers 1-7, 

PI and P2. Electronic reservoirs are denoted S, D (both blue) and H (red), (b) Stability diagram 
of the QD-system showing the conductance of QD2. The characteristic honeycombs are indicated 
with red lines. QD occupation numbers are denoted with N, M. AVc indicates the capacitive 
coupling energy, (c) Gurrent signal A/d in reservoir D with VsD ~ —30 fiV for Th ~ Ts,d + 100 mK. 
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FIG. 2. (a) Color scale plot of the change of current A/d with Th in the region of the honeycomb 
vertex. Black, dashed lines indicate the boundaries of the stability regions denoted 1-4. Data 
taken for VsD = —100 //V. Top and side panel show single traces extracted for Fpi = —381 mV 
and Vp 2 = —521 mV, indicated by horizontal (green) and vertical (red) line in the color scale plot 
(b) Data taken for inverted bias voltage: VsD = 100 /rV. (c) Conductance stability diagram for 
the same gate voltage region. 
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FIG. 3. Schematic energy diagram of the QD-system showing the alignment of chemical potentials 
for stability regions (a) [N + 1, M), section 1 and (b) {N +1, M +1), section 2. Each configuration 
is shown for low temperature (left) and high temperature (right) in reservoir H and with finite I^d- 
Red arrows indicate enhanced occupation fluctuations. 
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FIG. 4. Model calculations for Id and A/d- The following parameters were used: Ts^d = T^ase = 
230 mK, fig = —/td = —50 /reV, = 0 and Eq = 90 ^eV. (a) Id for Th = 230 mK (left) 
and Th = 330 mK (right), (b) Subtraction of the figures given in (a) yields A/d = /d(Th = 
330 mK) — /d(Th = 230 mK). 
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